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been  engaged  in  analysis  for  a ^ ^ ^ C°rporatl°n  has 
criteria  employing  concepts  oTboulT10^^  ^ hydr°dynamic  d^ign 
ing,  suction,  and  heating.  ^ C°ntro1*  including  shaj 

F°r  3 body  moving  at  lower  speeds  iam- 
before  considered  may  arise.  Th  . ’ fl°W  phenomena  never 

interaction  between  the  buoyancy^/  6001116113  ^ aSS°Ciated  with  the 

differences  and  a basic  bonndary  lays"  f"  ^ S'"a11  te,"Perat;ure 

0f  the  b°dy-  »*■«  analytic  methods  thl  8y"'MitriCal  ab°«  am, 
bleed  effects  of  buoyancy  and  forced  'flow^vTJV0^1""8 
a Preliminary  step  i„  the  delineation  f Oflloder,  as 

configuration  that  is  otherwise  08088  fl°“  effects  ""  a flow 

This  report  should  be  u.efulTThT  d'’OUt  'he  °f  the  b°d>’- 
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t0  the  predl«‘o„  and  improvement  of  the  perfoT"  yl"S  n“ld 
vehicles.  Performance  of  underwater 
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* SUMMARY 

buo  Heat  transfer  and  Shear  streSs  may  be  significantly  affected  by 
conlTtio  °tTS  anVSSOClated  ‘K-  »tlo„s  ln  many  forced. 

r iL8A;::;j10::;  “ “hen  * — — 

8 neated*  semi-infmite  horizontal  cylinder  Tho 
flows  doZsllllZZ  llZZuyll  ^ ^ ^ ^ 

~ for  moderate-speed  77“  “ 

know  the  development  of  k important  to 

eveiopment  of  the  buoyancy-force  effecr  ru  . 

which  the  effects  of  f,  conditions  under 

°f  fr«  convection  may  be  ignored  ln  heat-transfer 

rf:h“r“  CalCUlatl°nS*  a"d  ‘ha  beyond  uhich  tbe  effects 

fr;:  ~ion  may  te  as  in™  - *•  *«c-™llon  7 

scribed  b reP°rt’  ‘ S°1Uti°"  la  *«-  *«  tbe  flo„  de- 

bed  above.  An  asymptotic  expansion  based  on  tbe  ratio  of  the 

;;rf  rr  “ ^ — - -e  Reynolds  number  indicates  the 

gion  .here  buoyancy  forces  can  be  treated  as  second  order.  Tbe  ordi- 

7Z  tir^T  eqUatl°nB'  redUCed  by  a strailarity  transformation 
he  three-dimensional  boundary-layer  equations,  are  solved  by 

numerical  i„te8rati„n  over  tbe  Prandtl  number  range  from  0.01  to  10 

the  r suUs  supply  a quMtltative  criterion  ^ dlstinguifihing  ' 

onditions  under  which  the  effects  of  frPP 
neglected.  “ convection  may  not  be 

lOO’/h6  !T ‘T'  "hen  aPPlled  “ the  aXlal  “««■  1 - -ter  of  a 

eff  ct  "V  1CylIn<‘er  "ltH  1_ft  radi"s-  indicates  that  free  convection 
fects  should  be  considered  for  axial  rH 
eter  or  r k distances  greater  than  a diam- 

than  fi  d7  6 SP6ed  18  l0W  (~  4 ft/S6C)  ^ 3t  dlsta—  greater 

ban  f1Ve  diameters  at  higher  speeds  (~  20  ft/sec).  Furthermore,  the 

" eff7  lnCreaS6S  “ the  ^ ^ <=**-«  increases, 

secondary  flow  profile  is  required  to  be  small  for  tbe 

point  “rr  “ aPPly’  “ al“ayS  6Xhiblta  an  ln£laa“»" 

ect  of  this  on  boundary-layer  stability  and  transition 

in  Potentially  significant  and  should  be  investigated  further. 
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SYMBOLS 


a = radius  of  the  cylinder 

g = gravitational  acceleration,  Eq.  (3) 

g = similarity  nondimen si onal  temperature,  Eq.  (7) 

G = similarity  nondimensional  temperature  independent  of  <j>,  Eqs. 
(9)  and  (10) 

Gr  = Grashof  number,  Eq.  (3) 

f = nondimensional  stream  function,  Eqs.  (5),  (6),  and  (7) 

F = similarity  nondimensional  stream  function,  Eqs.  (9)  and  (10) 
Nu  = Nusselt  number,  Eq.  (13)  and  (14) 

Pr  = Prandtl  number,  Eq.  (3) 
r = radial  coordinate 
Re  = Reynolds  number,  Eq.  (3) 

T = temperature 
u = axial  velocity 
v = circumferential  velocity 
w = radial  velocity 
x = axial  coordinate 
a = thermal  diffusivity 
8 = thermal  expansion  coefficient 
e = Gr/Re^,  Eq.  (3) 

H = r//2ic,  Blasius  similarity  variable,  Eq.  (5) 

0 = nondimensional  temperature,  Eq.  (3) 

V = kinematic  viscosity 

<p  = circumferential  coordinate 

T = shear  stress 

Subscripts 

0 = zeroth  order  solution 

1 = first  order  solution 


r 

"t  *1 
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w = surface 
00  = free  stream 
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S uperscr ipts 

— = dimensional  quantities 
' = derivative  with  respect  to  n 
x = local  quantities  at  location  x 
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I ■ INTRODUCTION 


Free-forced  convection  has  been  a subject  of  great  importance  in 
such  engineering  devices  as  solar  collectors  and  nuclear  reactors  and 
in  such  applications  as  boundary-layer  stabilization.  The  effect  of 
free  convection  on  two-dimensional  boundary  layers  has  been  extensively 
studied  because  of  its  simple  nature— for  example,  the  heated  (or 
cooled)  vertical  flat  plate  with  upflow  or  downflow.  The  relative 
importance  of  free  convection  in  forced  convection  depends  on  the 
ratio  Gr-/Re-  = BgxAT/u^  , which  can  be  small  for  a short  plate  and 
becomes  one  of  the  dominant  effects  for  a large  plate.  Acrivos  [1] 
established  the  ranges  of  different  regions  using  an  integral  approach. 
Similar  results  were  obtained  by  Sparrow  and  Gregg  [2].  A more  com- 
plex two-dimensional  problem  of  co-flow  or  counter-flow  was  investi- 
gated by  Gersten  and  KBrner  [3].  They  looked  at  stagnation  flow  with 
suction  or  blowing  under  mixed  free-forced  convection  conditions.  The 
buoyancy  forces  were  shown  to  act  as  an  adverse  pressure  gradient  to 
counter-flow  and  to  lead  to  an  earlier  separation. 

Cross-flow  buoyancy  forces  may  be  generated  by  a heated  horizontal 
surface.  Under  conditions  which  do  not  lead  to  thermal  instabilities, 
analyses  have  been  reported  by  Mori  [4],  Sparrow  and  Minkowycz  [5], 
Hauptmann  [6],  and  more  recently  by  Redekopp  and  Charwat  [7].  It  was 
found  that  the  buoyancy  effect  had  to  be  considered  if  Gr_/Re_"^^  was 
greater  than  0.01.  The  heated-from-below  situation  causes  a more 
dramatic  effect  than  the  heated-from-above  situation.  This  is  because 
the  heated-from-below  case  results  in  an  inflection  point  in  the  ve- 
locity profile  and  magnifies  separation  as  well  as  initiating  thermal 
instabilities . 

"*  A hot  horizontal  cylinder  in  cross-flow  has  been  the  subject  of 
several  investigations.  The  most  recent  work  of  Oosthuizen  and  Hart 
[8]  presents  a criterion  for  consideration  of  buoyancy  forces  which 
depends  on  the  direction  of  the  flow  relative  to  the  vertical.  A 
vertical  heated  cylinder  in  upflow  or  downflow  has  also  been  investi- 
gated [9]. 
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An  exact  solution  of  the  three-dimensional  laminar  boundary  layer 
over  a heated  surface  when  gravity-driven  buoyancy  effects  introduce 
significant  cross  flow  into  an  otherwise  axially  symmetric  flow  geom- 
etry is  presented  in  this  report.  The  effect  of  surface  heating  on 
laminar  boundary-layer  stability  and  transition  has  been  analyzed  by 
Wazzan,  Okamura  and  Smith  [10],  and  the  effect  on  separation  by  Aroesty 
and  Berger  [11].  The  buoyancy  cross-flow  effects  have  not,  however, 
been  included  in  such  studies. 

The  physical  model  chosen  for  study  is  a semi-infinite  cylinder 
of  radius  a,  which  is  aligned  with  its  axis  parallel  to  a uniform 
flow  and  normal  to  the  direction  of  gravity.  The  uniform  flow  is 
assumed  to  have  a velocity  um  and  temperature  Tm.  The  surface  of  the 
cylinder  is  heated  to  a constant  temperature  Tw  (Ttf  > Tj . For  most 
external  flows,  the  buoyancy  force  can  be  neglected  in  a small  pure 
forced  convection  region  of  sice  0 (u„v1/3/ [Bg<T,  - TJ 1 ) dovmstre.™ 

of  the  leading  edge.  Beyond  that  region  the  effect  of  buoyancy  cross- 
flow  increases  as  the  fluid  flows  downstream.  There  is,  however,  a 
region  where  it  is  still  small  and  can  be  treated  as  second  order. 
Further  downstream,  a distance  of  order  aRe/Gr1/2,  the  initially  small 
buoyancy  effect  becomes  as  important  as  th.  forced  convection  effects 
and  can  no  longer  be  treated  as  a second  order  effect.  The  solution 
of  the  intermediate  region  where  the  buoyancy  effect  is  second  order 
is  presented  in  this  paper.  In  this  region,  the  fluid  which  is  en- 
trained by  the  basic  forced  laminar  boundary- layer  flow  is  given  a 
cross-flow  component  by  the  buoyancy  effect  of  the  heated  walls.  .Be- 
cause of  this  entrainment,  the  secondary  flow  velocity  grows  linearly 
in  the  downstream  direction.  This  is  a different  situation  from  the 
usual  pure  free  convection  flow,  where  all  boundary-layer  entrainment 
is  through  buoyancy,  and  where  the  axial  entrainment  is 


zero . 
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II.  ANALYSIS 


The  nondimensional  Boussinesq  boundary-layer  equations  in  cylindri 
cal  coordinates,  as  shown  in  Fig.  1,  are 


Fig.  1 — Physical  model  and  coordinates 


after  neglecting  smaller  order  terms  under  the  condition  that 


Gr/R«3/2  > 1 


(2) 


The  nondimensional  variables  used  in  scaling  Eqs.  (1)  are 


u = — , v = 


w 


w /Re 


(the  velocities)  ; 


e = 


T - T 

oo 

T - T~ 
w 00 


(the  temperature)  ; 


x = - , r = 


(r  - a)  /Re 


(the  coordinates)  ; 


Re  = 


u a 

oo 


V 


(the  Reynolds  number)  ; 


(3) 


6gaJ (T  - T ) 


Gr  = 


w 


(the  Grashof  number)  ; 


v 


Pr  = v/a  (the  Prandtl  number)  , and 


e = Gr/Re  , 


where  a is  the  radius  of  the  cylinder,  V is  the  kinematic  viscosity, 
u and  T are  the  free  stream  velocity  and  temperature,  respectively, 

OO  OO 

Tw  is  the  wall  temperature,  a is  the  thermal  diffusivity,  g is  the 
gravitational  acceleration,  and  B is  the  thermal  expansional  coeffi- 
cient, which  is  related  to  density  by  p = p^ll  - B(T  - T^) ] . 

As  mentioned  above,  for  the  region  very  close  to  the  leading  edge 
of  the  cylinder,  the  buoyancy  effect  is  negligible.  The  magnitude  of 
this  region,  defined  by  Eq.  (2),  is  of  the  order 

J ~u<aV1/3/[Bg(T  - TJ]2/3  . 


Downstream  of  the  leading-edge  region,  the  soxution  of  equations  .,n  (1) 
can  be  expanded  into  a series  of  e,  if  s is  small,  so  that 


u - Uq  + £u^  + 


(4-a) 


v = ev^  + 


(4-b) 


w = wQ  + + . . . 


(4-c) 


9 = 0Q  + e01  + ... 


(4-d) 


Substitution  of  the  expansion  given  by  Eqs.  (4)  into  Eqs.  (1), 
and  the  collection  of  terms  of  equal  order  will  result  in  the  perturba- 
tion equations.  The  perturbation  equations  of  lowest  order  are 


(e° ) : f"'+ff*  = n 

V ' 0 0 0 u ’ 


(5-a) 


V + PrfoV  - 0 - 


(5-b) 


where  the  prime  denotes  a derivative  with  respect  to  n,  with  r)  = r//2x 
being  the 
fined  by 


being  the  Blasius  similarity  variable.  The  stream  function  fQ  is  de- 


U0  f0  ^ ’ 


W°  ~ /2x  (nf°  ” f°} 


(5-c) 


(5-d) 


The  solution  of  Eq . (5-a)  is  the  Blasius  solution.  Equation  (5-b)  is 
the  forced-convection  energy  equation,  whose  solution  was  first  given 
by  Pohlhausen  in  1921. 

The  second  order  perturbation  equations  are 


(^):  f/ + f0V  - *f0'  V+5Vfl+  Vaf  = °’ 


(b 


V"  + foV  ~ V f2  + 4,  • eQ  = 0 , 


(b 


3f . 


+ Pr  ( fQg'  - 4f0'  g + 50Q'  f1  + 0Q'  ■—)  = 0 . (b 


The  stream  functions  f.^  and  f 2>  and  the  temperature  g are  defined  by 


ux  = (2x) 2 f1/(n,4>)  , 


v = (2x)  in‘  (r},<|>)  , 

1 


“l-  <2-)3/2(nf1'  -5^-^). 


q1  = (2x)  g(n,4»)  • 


The  boundary  conditions  for  Eqs.  (6)  are 

f^O.t)))  = f1/(o,<j>)  = f^  (°°,<j>)  = 0 , 

9f  9 

f2'(o,«j>)  = 3^  (0,4))  = f2'(°°,<j>)  = f2/(n,0)  = o , 


> < 


g(0)  = g(°°)  = 0 . J 

Equations  (6)  are  separable  in  terms  of  n and  (j)  witn  dependent  vari- 
ables of  the  form 


f1(n,4>)  - (n)  cos  (j)  , 

f2(n,<J0  = F2(n)  sin  <|>  , 


( 


g(n,4>)  = G(n)  cos  <t>  . 


Equations  (6)  can  be  further  reduced  to  the  following  ordinary 
differential  equations: 


V+foV-  V V • 


0 F1  + 5f0*  Fx  - -ffl'  F2  , 


(10-a) 


V+foV-  V V = -9o 


(10-b) 


G + (f0G  4fQ/  G)  = ~60/  (5F!  + V . 


(10-c) 


The  associated  boundary  conditions  for 


Eqs.  (10) 


Fx(0)  = ^'(0)  = Fl'(«)  = o , 
f2(0)  = F2'(0)  = F '(co)  = o , 


G(0)  = G (°°)  = 0 . 


olutions  of  Eqs.  (10)  are  integrated  numerically  for  Pr  = 0 01 
1,  8.  these  results  will  be  utilised  below  in  heat  transfer  and  shear- 
calculations.  It  is  of  interest  to  point  out  that  Eqs.  (7) 
indicate  that  the  effect  of  cross-flow  becomes  stronger  as  the  fluid 
moves  downstream.  Combining  Eqs.  (4),  (5-c),  (5-d),  and  (7)  gives 


u - V + £(2x)2  Fj'  cos  t + 


V = e(2x)  F ' sin  <J>  + ...  , 


W = /Tx  (nF°'  " Fq)  + E(2X)3/2  (nV  " 5Fi  + F2)  cos  <p  + ...  - 


® 90  + e(2x)^  G cos  <f>  + . . . . 


The  above  expansion  breaks  down  at  a distance  of  0(a  • Re/Gr  ), 
corresponding  physically  to  the  point  beyond  which  the  buoyancy  cross- 
flow  effect,  initially  small,  becomes  as  important  as  the  forced- 
convection  effects.  Solutions  for  the  region  beyond  0(a  • Re/Gr “)  will 
not  be  considered  in  this  report. 
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III.  RESULTS  AND  DISCUSSION 


Numerical  values  of  the  functions  F^,  , F^,  F * , and  G are 

presented  in  Figs.  2,  3,  and  4 for  Pr  = 0.01,  1,  10.  Equations  (12) 
show  that  the  buoyancy  force  stabilizes  the  flow  and  accelerates  its 
speed  over  the  lower  half  of  the  cyl^rder  (-tt/2  < < rr/2)  . Over  the 

upper  half  of  the  cylinder  (tt/2  < (J>  < 3tt/2),  the  axial  flow  is  de- 
celerated and  may  be  destabilized.  The  cross  flow  is  accelerated  by 
the  buoyancy  force  from  the  lower  stagnation  point  (<j>  = 0)  to  attain 
its  maximum  value  at  <f>  = tt/2,  and  then  decelerated  to  its  upper  stag- 
nation point  (<p  = n)  . It  is  interesting  to  note  that  the  solution  (12) 
represents  the  case  of  a vertical  flat  plate  whe  n <p  = tt/2  for  the  region 
near  the  leading  edge  where  the  cross  flow  is  small.  For  a vertical 
flat  plate,  the  cross  flow  is  decoupled  from  the  axial  flow;  in  other 


V (for  FV; 

4 3 2 10 


Fig.  2 — F|  and  F i 1 functions 
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words,  the  cross  flow  does  not  affect  the  axial  flow  and  the  tempera- 
ture distribution.  The  magnitude  of  the  circumferential  velocity,  v, 
is  being  increased  as  a linear  function  of  the  axial  coordinate,  and 
eventually  becomes  one  of  the  dominant  velocity  components  when  the 
fluid  flows  far  enough  downstream. 

The  local  Nusselt  number  with  respect  to  a and  (T_  - Tm)  can  be 
derived  from  Eq.  (12),  which  is 


w 


Nux  = "V?  [V  <0)  + e(2x)2  G'(°>  ’ cos  $ + ...]  . (13) 


The  first  term  of  Eq.  (13)  represents  the  Nusselt  number  for  a pure 
forced  convection.  The  relative  importance  of  heat  transfer  due  to 
cross  flow  can  be  indicated  by 


Nux/(Nux)fc  = 1 + e(2x)2  G1#(0)/60/  • cos  <J>  + ...  , 


(14) 


where  fc  denotes  the  forced  convection.  Equation  (14)  shows  that  the 
cross  flow  enhances  the  heat  transfer  on  the  lower  half  of  the  cylinder 


and  degrades  it  on  the  upper  half  of  the  cylinder.  Values  of  G7 (O)/0^/(O) 


are  given  in  the  following  table  for  Pr  = 0.01,  1,  10.  Equation  (14) 


Pr 

G'  (0)  /Oq'  (0) 

F1"(0)/f0//(0) 

f2"(o) 

0.01 

0.08638 

0.19509 

0.92627 

1 

0.09959 

0.09959 

0.61904 

10 

0.06576 

0.04458 

0.39761 

has  been  plotted  on  Fig.  5 for  <})  = 0 and  IT.  The  curves  clearly  show 
the  developing  cross-flow  effect  along  the  x-axis  and  show  where  the 
originally  small  cross-flow  effect  grows  and  will  eventually  overwhelm 
the  pure  forced-convection  one,  beyond  which  the  expansions  (12)  are 
no  longer  valid. 


. n 


i ! 


i 
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Fig.  5 — Heat  flux  diitribution 


The  local  shear  stress  at  the  cylinder  surface  can  be  computed 
from  the  equation 


and 


= y 


3v\ 

dr  n 
'r=0 


Introducing  the  series  expansions  (12),  the  relative  importance  of  the 
cross-flow  effect  on  the  axial  shear  stress  can  be  found  by 

T /(T  K = 1 + e(2x) 2 F/(0)/f  *(0)  • cos  <f>  + ...  . (15) 

rx  rx  fc  1 U 

The  circumferential  shear  stress  can  be  shown  to  be  proportional  to 

T . ~ e/2x  F /(O)  • sin  <f>  . (16) 

r<p  2 

Values  of  F1#(0) /fQ"(0)  and  F 2"  also  are  given  in  the  table  above.  The 
ratio  given  in  Eq.  (15)  has  been  plotted  in  Fig.  6 for  <}>  = 0 and  <f>  = IT. 
As  can  be  seen,  for  the  4>  = 71  curves  the  ratio  tends  to  approach  zero, 

r\ 

for  large  enough  ex  , which  indicates  that  separation  is  imminent.  For 
moderate  free— stream  velocities,  this  can  occur  within  several  radii  of 


the  leading  edge.  Further,  the  unevenly  distributed  axial  shear  stress, 

Eq.  (15),  can  induce  a pitch  moment  and  could  cause  an  oscillation  of 
the  moving  cylinder. 

Equations  (14)  and  (15)  indicate  that  the  cross-flow  effect  on 
heat  transfer  and  shear  stress  for  a heated  horizontal  cylinder  grows 
rapidly  when  the  fluid  flows  downstream,  proportionally  to  x2.  This 
means  that  an  initially  small  cross-flow  effect,  which  may  be  neglected 
in  the  region  close  to  the  leading  edge  of  the  cylinder,  cannot  be  ig- 
nored for  a heated,  long  slender  body. 

The  analysis  is  strictly  applicable  to  the  flow  over  the  outer 
surface  of  a heated  hollow  cylinder.  It  is  also  relevant  to  the  flow 
over  a slender  body  of  revolution,  provided  that  the  pressure  gradient 
associated  with  nose  shape  is  negligible.  The  effect  of  pressure  grad- 
ient and  nose  shape  will  be  reported  later. 

The  parameter  e has  been  evaluated  for  water,  subject  to  a 30°  wall 
overheat  and  a cylinder  radius  of  1 ft.  When  is  3 ft/sec  > e ~ .01 

and  when  u^  is  30  ft/sec,  e ~ .0001.  The  magnitude  of  the  cross  flow 

2 

is  proportional  to  ex  ; this  suggests  that  the  present  analysis  is 
valid  when  x < 10  ft,  at  low  speeds.  Within  this  region,  the  effects 
on  laminar  heat  transfer  and  shear  stress  are  not  overwhelming;  however. 
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Using  analytical  methods,  this  report  considers  the 
effect  of  forces  acting  counter  to  laminar  flow 
around  heated  bodies  moving  through  water  at  low 
speeds.  A cross  flow  is  induced  when  a uniform 
horizontal  stream  passes  along  a heated  semi- 
infinite  horizontal  cylinder.  The  cross-flow 
effects  on  heat  transfer  and  shear  stress  grow  as 
the  fluid  flows  downstream  and  eventually  becomes 
one  of  the  dominant  flow  mechanisms.  When  the 
analysis  applied  to  specifically  a cylinder  of 
one  foot  radius,  heated  to  100  degrees  F,  at 
speeds  of  about  four  feet  per  second,  free  con- 
vection effects  should  be  considered  for  distances 
as  great  as  four  feet  from  the  leading  edge;  at 
about  20  feet  per  second,  out  to  ten  feet  from 
the  leading  edge.^  (BK) 


